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Abstract — Consider a two dimensional domain 5 C 
containing two sets of nodes from two statistically independent 
uniform Poisson point processes with constant densities pl and 
Pnl- The first point process identifies the distribution of a set 
of nodes having information about their positions, hereafter 
denoted as L-nodes (Localized-nodes), while the other is used 
to model the spatial distribution of nodes which need to localize 
themselves, hereafter denoted as NL-nodes (Not Localized-nodes). 
For simplicity, both kind of nodes are equipped with the same 
kind of transceiver, and communicate over a channel affected by 
shadow fading. 

As a first goal, we derive the probability that a randomly 
chosen NL-node over S gets localized as a function of a variety 
of parameters. Then, we derive the probability that the whole 
network of NL-nodes over S gets localized. 

As with many other random graph properties, the localization 
probability is a monotone graph property showing thresholds. 
We derive both finite (when the number of nodes in the bounded 
domain is finite and does not grow) and asymptotic thresholds 
for the localization probability. 

In connection with the asymptotic thresholds, we show the 
presence of asymptotic thresholds on the network localization 
probability in two different scenarios. The first refers to dense 
networks, which arise when the domain 5* is bounded and the 
densities of the two kinds of nodes tend to grow unboundedly. 
The second kind of thresholds manifest themselves when the 
considered domain increases but the number of nodes grow in 
such a way that the L-node density remains constant throughout 
the investigated domain. In this scenario, what matters is the 
minimum value of the maximum transmission range averaged 
over the fading process, denoted as dmax, above which the 
network of NL-nodes almost surely gets asymptotically localized. 

Index Terms — Ad-hoc network, connectivity, GPS, LBS, local- 
ization, location based services, positioning, probabilistic method, 
random arrays, sensor networks. 

I. Introduction and Literature overview 

This paper deals with a network composed of two sets of 
nodes randomly distributed over a two dimensional domain 
5 C SR^ following two statistically independent Poisson point 
processes with intensities and pnl- The first process is 
associated with the nodes that have a-priori knowledge about 
their position (these are the so called L-nodes), while the other 
point process is associated with the nodes that are trying to 
localize themselves (these are the so called non-localized or 
NL-nodes). In particular, the paper focuses on the connection 
between some system level parameters and the node local- 
ization probability in a Poisson distributed configuration of 
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nodes, which are at the basis of topological network control. 
We do not propose any new or modified localization method. 
As it will become clear later, the primary assumptions in 
our analysis are: a) nodes are Poisson distributed over a 
bounded circular domain contained in 'St^ and b) each node has 
an average typically circular footprint representing its radio 
coverage. Hence, while we focus on a particular example 
involving range measurements using Received Signal Strength 
(RSS), the analysis can be applied to other range measurement 
methods as well. Notice that Poisson point processes are useful 
for modelling scenarios in which the deployment area, the 
number of deployed nodes, or both, are not a-priori known. 
The Poisson model is in fact a good approximation of a 
binomial random variable when the number of deployed nodes 
over a bounded domain is high while the node density is 
constant across the whole region of interest [1]. Nevertheless, 
the Poisson approximation leads in many cases of interest to 
a mathematically tractable problem. 

This general framework can be recognized in many practical 
scenarios. A possible example is a Distributed Sensor Net- 
work (DSN), in which one may be interested in distributed 
power efficient algorithms to derive localization information 
in a randomly distributed collection of severely energy and 
computation power limited nodes. A second example may 
be that of a wireless network, in which the various network 
elements may communicate between themselves (in the case 
of wireless networks allowing peer-to-peer communication) or 
with a subset of nodes whose positions are known (this is 
the case of classic cellular networks and WLANs, whereby 
every node must communicate with at least one base-station 
or access point). With this scenario in mind, let us provide 
a brief overview of the localization methods that have been 
proposed in the literature. 

Given the great difference between the communication 
and computation capability of the nodes, as exemplified by 
the DSN and WLANs, algorithms developed for localization 
should be tailored to the particular scenario at hand [2], [3]. 

Practical localization algorithms can be classified in at 
least two ways: centralized or distributed [2] and range- 
free or based on ranging techniques [4]. The most common 
techniques are based on measured range, whereby the location 
of nodes are estimated through some standard methods such as 
triangulation. Cramer-Rao Bounds (CRBs) on the variance of 
any unbiased estimate based on the above ranging techniques 
are readily available and provide a benchmark for assessing the 
performance of any given algorithm [5], although we should 
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note that the derivation of the CRB itself relies on a proba- 
biUstic model (often assumed to be Gaussian), that describes 
the connection between the parameter to be estimated and the 
raw observations. 

In range-free localization, connectivity between nodes is a 
binary event: either two nodes are within communication range 
of each other or they are not [6]. For simplicity, we may view 
this event as obtained from hard quantization of, for instance, 
a RSS random variable. If RSS is above a certain detection 
threshold, the nodes can communicate, otherwise they cannot. 
Of course, the nature of path loss and the terrain characteristics 
influence both the coverage radius and the deviation of the 
coverage zone from the ideal circular geometry. In a typical 
scenario there may be multipath. Multiple Access Interference 
(MAI) and Non Line Of Sight (NLOS) propagation conditions 
[2]. Various range free algorithms have been proposed in the 
literature including the centroid algorithm [7], the DV-HOP 
algorithm [8], the Amorphous positioning algorithm [9], APIT 
[10], and ROCRSSI [4]. 

A review of various localization techniques proposed in the 
literature may be found in [11]. In [12], the authors propose 
an approach based on connectivity information for deriving 
the locations of nodes in a network. In [13], the authors 
present some work in the field of source localization in sensor 
networks. 

A topic somewhat related to the problem dealt with in 
this paper is network connectivity. This topic has received 
much attention recently [14], [15]. Given n homogeneous 
nodes independently and uniformly distributed over a region 
5 C SR^, a network is said to be connected if there exists a 
communication link between every pair of nodes in S. Early 
work on this topic can be found in [16], [17], [18]. 

In [16], the authors investigated the percolation of broadcast 
information in a multihop one-dimensional radio network 
modeled by a Poisson spatial process. In [17], [18], the authors 
investigated the connectivity of two and one dimensional 
networks respectively, as a function of the transmission range 
of the nodes involved in the network. 

The seminal work [19] by Gupta and Kumar demonstrated 
that a network constituted by n i.i.d. randomly distributed 
sensors over a disk of area S, is asymptotically (i.e., for 
n oo) almost surely connected if the transmission range 
between nodes is chosen as 

r{n) = S ■ (log(n) + 7(n))/(7rn) 

provided that j{n) ^ oo as n oo. A more careful look 
at the asymptotic expression for r{n) above would reveal a 
resemblance to a known result on random graph theory [20] 
which states that given a set of n nodes, the random graph 
formed by adding an edge between any couple of nodes with 
probability p{n) will become connected almost surely if 

p{n) = (log(n) + j{n))/n 

as n ^ oo, provided that 7(n) — > oo as n oo. 

In [21] Xue and Kumar demonstrated that in a random 
network of n homogeneous nodes, the number of neighbors 
of a randomly chosen node required for the network to be 
asymptotically connected is 0(log(n)) as n ^ oo. Such 
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Fig. 1 . Pictorial representation of a bipartite network witli an average number 
|5|pl of L-nodes and \S\pnl of NL-nodes over a bounded domain S with 
size |5| = ttR^. 

results have been extended to 3 -dimensional networks in [22]. 
Other works focusing on the connectivity of random networks 
over bounded domains may be found in [23] -[25]. Finally, 
paper [26] studies the connectivity of multihop radio networks 
in log-normal shadow fading environment by looking at the 
probability that a randomly chosen node is asymptotically 
isolated. 

The rest of the paper is organized as follows. In Sec- 
tion |II] we formulate the problem at hand, present the basic 
assumptions for the derivations that follow, and briefly recall 
the mathematical notation needed in connection with the 
evaluation of the asymptotic thresholds. Section |III] recalls 
the mathematical models adopted for the characterization of 
the transmission channel between the two kind of nodes. The 
localization probabilities are derived in Section|IV]for a variety 
of transmission parameters. Section |V] investigates the pres- 
ence of finite thresholds above which the derived localization 
probabilities manifest large variations. This analysis is then 
extended in Section [Vll taking into account the behavior of the 
localization probabilities for unboundedly increasing values of 
the number of deployed nodes. Finally, Section [VTIl is devoted 
to conclusions. 

II. Problem Formulation and Assumptions 

Consider a circular domain C of radius R and 
area |5| = ttR^ where sensors are deployed following two 
statistically independent two dimensional Poisson point pro- 
cesses with uniform densities and pnl, respectively. For 
simplicity, both L and NL-nodes are assumed to employ the 
same kind of receiver and communicate in a scenario whereby 
the transmission channel is affected by shadow fading with 
variance a^. Two nodes can communicate if the received 
power is above a prespecified threshold Pw,th, which is 
a network parameter with respect to which the results are 
derived. 

L-nodes have localization information relative to some co- 
ordinate frame. Notice that how this localization is established 
is irrelevant to our problem formulation. 

On the other hand, NL-nodes need to localize themselves. 
Since we have two kinds of nodes, the connection model 
between them can be specified as a bipartite random network, 
denoted by Gl,jvl(PLi Patl)- A pictorial representation of a 
bipartite graph is shown in Fig. [T] whereby an edge between 
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the j-th NL-node and the i-th L-node is used to identify a 
communication Unk between the underUned nodes. Owing to 
the constant densities and pnl, the average number of L 
and NL-nodes over S is, respectively, ■ \S\ and p^L ■ \S\. 

The localization problem is two dimensional and three 
distance measurements relative to nodes with known positions 
are sufficient to solve for the {X, Y) coordinates of the NL- 
node unambiguously. 

A. Notations 

Throughout the paper we assume the following notations 
[27]. 

• x{n) = 0{y{n)) if there exists a suitable constant c such 
that x{n) < cy{n) for any n > Uq- Notation x{n) = 0{1) 
is used to signify that x{n) is a bounded sequence. 

• x(n) — o{y{n)) if 

lim [ I = 
n-+oo y[n) 

• x{n) ^ y{n), i.e., x and y are asymptotically equivalent, 
if and only if 

,. ^in) , 
lim ——^ = 1 

n— >oo y[n) 

It is a matter of fact that the previous condition can also 
be represented as follows: 

x{n) = y(n) + o{y{n)) = y{n){l + o(l)) 

• An event El which depends on the integer- valued vari- 
able N is said to be asymptotically almost sure (a.a.s), 
or to occur with high probability (w.h.p.), if 

lim P(El) = 1 

III. Random Graph Models for Wireless Networks 
OF Randomly distributed Nodes 

Connections between the two classes of nodes depend on 
the considered channel model. Basically, three basic models 
have been extensively adopted in the literature for wireless 
networks analysis, namely random geometric graphs [28], 
path-loss channel model [29], and path-loss geometric model 
with shadowing [29], [15], [26]. 

A. Random Geometric Graphs 

A random geometric graph suitable for the problem at hand, 
is defined as follows. Let 



xfi, xf^-z) identify the geometric 



position of the j-th NL-node, X^^, with j = 1, . . . , pnl\SI 
and let D = || • |j be some suitable nor on 3?^. In a random 
geometric graph, X^'" is connected to a L-node with 
i — 1, . . . , Pl\S\ over the domain S by an undirected edge if 
D — Xj^ll < r, whereby r is some positive predefined 

parameter 

This is a reasonable assumption in practice. In fact, usually 
receivers have strict signal-to-noise (SNR) requirements such 
that if the SNR is above a predefined threshold, i.e., if the dis- 
tance between the nodes is below a given value, then reliable 
communication between the nodes is possible; otherwise, no 
communication is allowed. 

'a thoroughly employed norm is the Euclidean norm. 



B. Path-loss Geometric Random Graph, Without Shadowing 

A somewhat better model accounting for practical commu- 
nication receivers is the so-called path-loss geometric random 
graph. 

Let us assume that the j-th NL-node can communicate with 
the z-th L-node if the power received by the i-th L-node is 
greater or equal to a certain threshold Pw.th- The coverage 
area of the j-th NL-node comprises the L-nodes where the 
received power from NL-node j is greater than or equal to 
Pw,th- A NL-node can only communicate directly with L- 
nodes that fall inside its coverage area. With this setup, we 
can model the presence of a communication link between the 
j-th NL-node and the i-th L-node with a random variable Ij ^ 
as shown in Fig.[T] Ij i is a discrete random variable assuming 
two possible values with probabilities Pji and 1 — Pji, i.e. 
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Based on the observations above, the probability Pji ~ 
P{Ij,i = 1) is equal to the probability that the power received 
by the i-th L-node is greater or equal to the power threshold 

Pw,th- 

Let us consider the power P{dj_i) received by the i-th L- 
node at a distance dj,i from the j-th NL-node [29]: 

{4n)^d]>^ 

whereby, Pt is the transmitted power, Gt is the transmitter 
antenna gain, G,- is the receiver antenna gain, rip is the path- 
loss exponent, and A = c// is the wavelength. Notice that this 



equation is not valid for dj 



0. 



The path-loss in dB PL-[dB] can be expressed as: 



PL[dB] = 10 log 



10 



Pt 



Pidj,r) 



= -10 log 



GtGr^'^ 



10 



(47r)2d" 



(2) 

Since this equation is not valid at dj^i = 0, usually it is 
specified with respect to a reference distance do. In other 
words, the received power P{dj,i) at a distance dj,i from 
the transmitter is given with respect to a reference power 
Po received at a distance dp, usually assumed equal to 1 
meter [29]. Such a value may be measured in a reference 
radio environment by averaging the received power at a 
given distance close to the transmitter Doing so, the equation 
specifying the received power P{dj i) is then expressed with 
respect to Po'. 



P{dj^i) = Po ■ 



do 



Po 



\/dj,i > do (3) 



, dj,i J \ do 

whereby Po is the signal power at a reference distance do 
normalized to one for simplicity, and Up is the path loss 
exponent. In a similar fashion, if we consider the receiver 
threshold power Pw.th, and define dmax as the distance 
between the transmitter and the receiver at which the received 
power P{dj.i) equals Pw,th, we can write: 



P{dj,i) — PiLKth 



dj,: 



— Pw.th 



di 



(4) 
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With this setup, the probabiHty Pji = P{Ij.i 
connection between a NL-node and a L-node can be evaluated 



1) of a Unk Upon setting a ~ 
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as: 



1, < djA < d, 
0, d, 



< R 



< dj^i < R 



(5) 



whereby R is the radius of the area on which the network is 
estabhshed. Notice that any distance must be smaller than R, 
and that in this model the radio coverage of any node is a 
perfect circular area with radius d„iax- Any L-node falling 
in a circle of radius d,nax from the NL-node is assumed 
to communicate with the reference NL-node. In this respect, 
dmax is the coverage radius of any node, and takes on the 
same meaning as r in the geometric random graph model 
described in the previous section. The difference is that here 
dmax is related to typical transmission conditions, while r 
in the previous section is only interpreted as a geometric 
parameter. 

The only parameter of interest in this model is the maximum 
distance dmax- Simulation results can be given with respect 
to the normalized distance in order to highlight the 

dependence of the results from the ratio between the coverage 
radius of any node and of the overall deployment area. 

C. Wireless Channel Model: Path-loss Geometric Random 
Graph with Shadowing 

Practical measurements of the signal power level received 
at a certain distance from a transmitter often indicate that 
the path-loss in (|2]i follows a log-normal distribution [29]. 
From one easily evaluates: 
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Let us consider the normalized variables P{dj^i) and dj_i, 
defined as 

The log-normal model is formalized as: 

lOlogio (PK-,0) = lOlogio [idJ,^)-'''-] +Xs 

whereby, Xs is a Gaussian-distributed shadowing random 
variable, i.e, Xg ~ N{fis, ) with /i^ = 0. With this setup, the 
probability that a NL-node and a L-node establish a wireless 
connection is: 

P(l01ogio {P{d,.,))>0) 

Notice that the underlying model becomes a path-loss geomet- 
ric random graph without shadowing upon setting ag — 0. 

By consMering P{dj,i)dB = lOlogi and yUd = 

lOlogiQ , it easily follows that: 

P {P{d,^i)dB > 0) = P (X, > -/id) 
The latter equation corresponds to: 
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can be rewritten as follows 



and rj = the previous equation 



P {P{d,^i)dB > 0) = 



1 - erf 



a 



ln(dj,i) 



(6) 



This is the probability of establishing a wireless link between 
a NL-node and a L-node given that their relative distance is 

dj,t- 

Let us focus on the bipartite graph of Fig. [H and assume 
that the j-th NL-node can communicate with the i-th L-node if 
the power received by the i-th L-node is greater than or equal 
to a certain threshold Pw,th- The coverage area of the j-th NL- 
node comprises the L-nodes where the power received from 
the j-th NL-node is greater than or equal to Pw,th- A NL-node 
can only communicate directly with L-nodes that fall inside its 
coverage area. However, with respect to the model described 
in the previous section, here there is a non-zero probability 
of a wireless communication between nodes that are far apart 
more than dmax due to the considered shadow fading model. 

With the setup above, we have: 
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whereby. 



Pth = 10 log 
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With this setup, we can model the presence of a communica- 
tion link between the j-th NL-node and the i-th L-node with a 
random variable ^ as shown in Fig. [T] The random variable 
Ij^i is a discrete random variable assuming two possible values 
with probabihties Pji and 1 — Pji like in ([T]), where 

Pj, - P (P{dj.,i)dB > 0) (9) 

as in (|6]). This is the most general model since when cTs — 
it becomes a path-loss geometric model. Moreover, upon 
assuming dmax — r, the geometric random graph described 
by Penrose [28] is obtained. 

IV. The Localization Probability 

The aim of this section is to derive the localization proba- 
bihty of the network of NL-nodes over the bounded domain 
S. The problem is solved by first determining the localization 
probability of a randomly chosen NL-node over S, and then 
upon identifying the localization probability of the set of NL- 
nodes falling within S based on justifiable assumptions. 

Owing to the definition of the Poisson point process de- 
scribing the NL-nodes distribution over S, the problem can 
be solved by evaluating the expected number Xnl.b. — 
E{d'^f^\R} of L-nodes seen by a NL-node within a circular 
area of radius R centered on the NL-node. Such a random 
variable is denoted as d^^. Resorting to ideas from percolation 
theory [30], the expected value of neighbors within a distance 
P of a generic NL-node can be evaluated as follows: 

^|P} = • / plP {P{r)dB > 0|r) rdrdcf) (10) 

Jo Jo 

whereby, p^, is the density of the point process related to the 
L-nodes, and P [P{r)dB > 0|r) is as defined in (|6]l with r — 
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Fig. 2. Behavior of the difference Xnl ^ Xnl,r as a function of the radius 
R of the considered domain S. All curves are related to =0.1 nodes/m"^. 
Other transmission parameters are as noted in the legend. 



The solution of (fTOl i. whose proof is reported in Appendix 
I, is: 



— R -TrpL— erf -In 
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(11) 



The expected number Ajvl = E{d'^f'"} of L-nodes seen by a 
NL-node over the entire can be evaluated as follows: 



lim / / pL-P{P{r)dB>0\r)rdrdcf> 
n^^Jo Jo ^^^^ 

The solution of (fT2l l. whose proof is given in Appendix I, is: 



A 



NL 



(13) 



Before proceeding further, notice that so long as i? ^ d„iax, 
the average number of L-nodes estimated by (fTTl i over S* C 3?^ 
coincides with the ones estimated by ( fT3] l over the whole two 
dimensional domain 3?^. This is clearly depicted in Fig. |2]as 
a function of the radius R of the considered domain S, for 
a variety of transmission parameters as noted in the legend. 
Actually, the less stringent condition R > 5 ■ dmax suffices 
to ensure Aatl « \nl.r- Owing to this observation, when 
not differently specified, in what follows we will consider the 
formula (fT3l i. 

The next line of pursuit consists in the definition of the 
localization probability of a randomly chosen NL-node within 
S. Since L-nodes are distributed as a Poisson point process, 
the number of L-nodes is a Poisson random variable with 
expected value \nl = E{d^^} in \f S = 3?^, or \nl,b. 
in ( fTTT l if 5 is a bounded domain of radius R contained in 
3?^. The event of interest, identified by El, is the event that 
a randomly chosen NL-node is within the transmission range 
of at least three L-nodes. Over 3?^, such a probability can 
be evaluated as the probability that the random variable 



Fig. 3. Behavior of the localization probability P{Ei^) as a function of the 
L-node density over R^. Other transmission parameters are as noted in 
the legend, while pjvz, = Pl- Simulated points are identitied by star-marked 
points over the respective theoretical curves. 



takes on values greater than or equal to 3: 

j 



P{El) = P(dr>3)^^ ^ ^ e-^l^" 

3=3 

" E{d^ 



= 1-E 

which can be rewritten as 



(14) 



PiEL) = l-e 
Using ( HjI ), it is straightforward to obtain: 



i+E{d:f^}+^^^^ 



+ ^7r2d4 



1 + PLTTdl^^^e^ + 



(15) 

The behavior of P [El] is displayed in Fig. [3] for the param- 
eters noted in the legend. 

Simulation results have been obtained as follows. We define 
a square domain C with size Rd x R^ and centered a 
circular domain S of area irR? in the middle of C. In order 
to simulate the entire domain 3?^, we assume Rd ^ R- 
Furthermore, we must have R ^ d^ax in the investigated 
scenario, say R > 10d„iax7 based on the considerations stated 
above. Then, we generate two statistically independent point 
processes distributed uniformly over C with densities and 
Pnl, respectively. Owing to the constant density of both point 
processes within C, the number of L-nodes falling in C is, on 
average, Ec — Pl - R% while the average number of L-nodes 
falhng in S is Er = pL ■ nR^ p^ ^ Er/ttR^. Upon 
substituting in Ec the following relation follows: 

R^ 

Ec = Er ■ 

For ensuring an appropriate number of L-nodes in S, say 
Er ~ 500, Ec nodes are uniformly distributed on the bigger 
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Fig. 4. Minimum L-node density over as a function of as (in dB) for 
assuring that on tlie average, each NL-node is able to establish a wireless link 
with at least three neighbors under the channel conditions exemplified by the 
parameters P^^ and rip. 



Fig. 5. Exp ected number AjvL = of L-nodes seen by a NL-node 

over 5R2 (see (HJ) as a function of the L-node density pj^ , for the transmission 
parameters noted in the legend. Simulated points are identified by star-marked 
points over the respective theoretical curves. 



domain C. The localization probability is then evaluated by 
dividing the number of localization events in the domain S 
by the number of randomly generated network realizations. In 
order to avoid border effects, NL-nodes close to the border 
of the domain S are allowed to communicate with L-nodes 
within an annulus of radius dmax from the circular domain S. 

Some observations from the results in Fig. [3] are in order. 
As expected, the node localization probability increases for 
increasing values of the density of the L-nodes. For fixed 
values of pL, the node localization probability increases for 
increasing values of the parameter Pth, which in turn depends 
on the maximum transmission range dmax- Moreover, note 
that for a given set of transmission parameters, the localization 
probability increases for increasing values of the variance of 
the shadow fading Cs. 

The analysis above is the starting point for finding theo- 
retical conditions assuring that the localization probability is 
above a certain threshold. Upon imposing E{d^^} > 3, one 
easily finds: 

2 

3 

PL > —n — e °' (16) 

which yields the minimum uniform L-node density over for 
assuring that on the average each NL-node is able to establish 
a wireless link with at least three neighbors under the channel 
conditions exemplified by the parameters as and rip. 

The behavior of (fT6] l as a function of the shadowing 
parameter as (in dB) is displayed in Fig.|4]for the transmission 
parameters noted in the legend. Notice that, as expected, 
shadowing tends to decrease the L-node density since farther 
nodes can communicate over longer distances. 

The behavior of the expected number Ajvl ~ E{d^^} of 
L-nodes seen by a NL-node over SR^ (see ( fT3b ) is displayed in 
Fig.|5]as a function of the L-node density p^ for a variety of 
transmission parameters, as summarized in the figure legend. 
Star-marked points denote simulated points. 

Next, consider the probability that the whole network of NL- 
nodes falling in the bounded domain S under investigation gets 



localized. Such an event occurs when all the single NL-nodes 
within S get localized. Let Nj^i^ be the number of NL-nodes 
falling within S. 

Consider P{El) in ( fTSl l. and define X{Xnl) as 

X(\nl) = 1 - P(£;l) = 6-"--'^™-'=^ [1+ 

7,2 2 ,,2 ] (17) 

With this setup, by virtue of the independence of the NL- 
nodes in S, the probability P/v (El) that all the network of 
NL-nodes deployed in S gets localized can be expressed as: 

Pn{El) = [I-X(A^l)]^"- (18) 

whereby, we have to interpret such a probability as conditioned 
on the number of NL-nodes falling in the domain S. On 
average, Nnl ~ Pnlt^R^ in the observation area S. 

V. Analysis of the Localization Probability and 
Thresholds, Finite Case 

Returning to our analysis where we assume the knowledge 
of the radio coverage area of a given NL-node, a common 
characteristic of many problems tackled using the probabilistic 
method is the existence of transition thresholds where the event 
of interest exhibits a large variation. Indeed, it is known that 
every monotone graph property in randomly generated graphs 
has a sharp transition threshold [31], [32]. Such thresholds 
are established in the asymptotic case, i.e., in the limit when 
the number of nodes in the random graph tends to infinity. 
Thresholds are very useful in practice for topology control of 
the network [14]. 

In what follows, we will first derive transition thresholds 
for the localization problem in finite regimes, i.e., when 
the numbers of both L and NL-nodes are finite within a 
bounded domain S as defined in the previous sections. In the 
second part, we will investigate the localization problem in 
the limiting cases of dense networks. Notice that our results 
hold even in the random geometric model by setting as — 0. 
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Fig. 6. Behavior of the localization probability P(i?i) as a function of 
the L-node density p^. Other transmission parameters are = 40 dB, 
(Tg = 4 dB, Up = 2, pML =0.1 NL-nodes/m^ and R = 100m. 



Fig. 7. Behavior of the localization probability PpfiEi^) as a function of 
the L-node density p^. Other transmission parameters are /^jj, = 40 dB, 
(Ts = 4 dB, Up = 2, pjvl = 0.1 NL-nodes/m^ and R = 100m. 



A. Thresholds for Single Node Localization Probability, Finite 
Case 

Since the localization probability P{El) in ( fTSb is a mono- 
tonically increasing function of its arguments embraced within 
Xnl^ the transition thresholds observable in the finite regime 
(especially for large values of p^) can be obtained by taking 
the second partial derivative of P {El) in ( fTSt with respect to 
the parameters of interest, such as pL and dmax, and setting 
the result to zero. 

Let S be the usual bounded circular domain of radius R in 
5R^. Let us analyze the thresholds of P (El) with respect to 

PL- Let 7i = T^df^ax^'^ ■ After some algebra, the first partial 
derivative with respect to can be expressed as 



d 

^P{El) 

OPL 



(19) 



Given that pL > 0, ( fT9] l is always greater than zero, showing 
a strictly increasing behavior of P (El) with respect to pL- 

The second partial derivative T{pl) = -^-^P{El) of 
P {El) with respect to pL is: 



T{pl) = e-^-P-^lpL 



1 ^1 



(20) 



The values of the threshold are the solutions of the equation 
T{pl) = 0, that is, 

2 



1 



71 



PL =0: 



TTdl, 



(21) 



Fig. |6] shows the behavior of the localization probability 
P{El) as a function of pL for the transmission setup noted 
in the figure caption. Moreover, in the figure we report the 
behavior of the second derivative T{pl) (normalized with 
respect to its maximum for depicting both curve on the same 
ordinate range) along with the threshold obtained by (ISTT i 
with the setup noted above. 

Let us analyze the thresholds of P {El) with respect to 
dmax, and for ease of notation, set d„iax = d„i and 72 = 
PLT^e^ . Following the same reasoning as applied for p^. 



after some algebra, one easily obtains the threshold for the 
localization probability with respect to the node transmission 
range 



(22) 



B. Thresholds for the Localization Probability of the Whole 
Network of NL-nodes, Finite Case 

Owing to the fact that Pn{El) < P{El) for a given 
transmission scenario, thresholds for the probability Pn{El) 
are expected to be higher than the ones obtained for P{El)- 

Let us start our analysis by deriving the thresholds of 

Pn {El) in JTSl l with respect to Pl- Let 71 — Trd^^^e^. 
After some algebra, the second partial derivative F{pl) = 
S-tPn {El) of Pn {El) with respect to pL is: 

Op^ 



F{pl) = WiNnl 1 - e 



-IiPl 



~^-P-{2pL~lipl)+^ 



l + llpL + ^ltpi) 

llp'ie-''-<^''i^(NNL~l) 



Nr. 



2 i^e-'<i''L{l+-,iPL+k'llpl) 

(23) 

The values of the threshold are the solutions of the equation 
F{pl) = 0. Noting that 



1 



llPL 



^ 2 2 

-,1iPl 



with PL > and 71 > 0, the only solutions are the roots of 
the non-linear equation: 



1 



llPL 



1iPl< 



-llPL 



{N. 



NL 



I) 



2 1 - e-TiPi (1 + 7iPl 



bfpl) 



(24) 



As a reference example, consider the transmission scenario 
investigated in the previous section, and summarized in the 
caption of Fig. [7] which shows the behavior of the localization 
probability Pn{El) as a function of p^. Also shown is the 
behavior of the second derivative F{pl) (normalized with 
respect to its maximum for depicting both curve on the same 
ordinate range). Note that the threshold for Pn{El) is about 



Fig. 8. Finite case thresholds logj^Q(p^) for Pn{El) as a function of fSt^ 
for a variety of parameters noted in the legend. Other transmission parameters 
common to all plots are Pml =0.1 NL-nodes/m^ and R = 100m. 
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Fig. 9. Behavior of the localization probability Pn{El) as a function of 
the maximum transmission range dmax- Other transmission parameters are 
/3(^ = 40 dB, (Ts = 4 dB, rip = 2, = pjvi = 0.1 NL-nodes/m^ and 
R = 100m. 



one order of magnitude greater than the threshold noted 
in dlB, relative to P{El). 

The behavior of the thresholds (obtained as the solutions 
of (l24l i) as a function of the parameter (3tfi for various values 
of the path-loss exponent Up and as is depicted in Fig. |8] 
From this figure, we observe the decreasing behavior of 
for increasing values of Pth, i-C- for increasing values of the 
maximum transmission range dmax noted in (|7]i. 

Let us analyze the thresholds of Pn {El) with respect to 
dmax, and for ease of notation, set dmax = d„i- Let 72 ~ 
Plttc'^ . After some algebra, the second partial derivative 
F{pl) ^ S^Pn (El) of Pn (El) with respect to dm is: 



F{d,n) = 1 - e-T^'*" (1 + 72^2 



2 12 "m ) 



Fig. 10. Finite case thresholds of the localization probability Pn{El) 
as a function of the L-node density pi^ for a variety of parameters noted in 
the legend. Other transmission parameters are pml =0.1 NL-nodes/m^ and 
R = 100m. 



The values of the threshold d* are the solutions of the 



equation F{dm) 



jlNNLdf,e-^^-< 



272^?, 



-1) 



(25) 



P+720 



0. Upon noting that 

^ „2j4 



> 1 + 72< + 7772^™ , Vd™ > 0, 72 > 



the only solutions are the roots of the non-linear equation: 

{Nnl - 1) 



272<, 



lldte 



= (26) 



l-e-^=''-(l + 72rf^ + hl^i) 

Fig. |9] shows the behavior of the network localization proba- 
bility Pn{El) as a function of dm for the transmission setup 
noted in the figure caption. The figure also shows the behavior 
of the second derivative F{dm) (normalized with respect to 
its maximum for depicting both curve on the same ordinate 
range) along with the threshold obtained by solving the 
non-linear equation ( |26] ) with the setup noted in the caption 

of Fig. m 

The behavior of the thresholds (obtained as the solutions 
of ( l26b ) as a function of the L-node density pL for rip = 4 and 
various values of as is depicted in Fig.[TOl From this figure, we 
observe the decreasing behavior of d\-^ for increasing values 
of PL- 

VI. Asymptotic Behavior of the Localization 
Probability and Thresholds 

In this section, we present results on the behavior of the 
localization probabilities of both single NL-node and the 
overall network of NL-nodes deployed over both bounded 
and unbounded domains in a transmitting scenario affected 
by shadow fading. 

The first result concerns dense networks, i.e., network of 
nodes whereby the node densities of both point processes 
deployed over a disk C with radius R ^ dmax, are 
allowed to grow unboundedly as a function of the number 
of nodes over S. As above, edge effects are neglected, and 
the hypothesis R ^ dmax allows us to employ the relation 
^NL.R ~ ^NL- Moreover, assume that the transmission range 



9 



is homogeneous and equal to dmax for both kinds of nodes. 
The next theorem investigates the behavior of the locaHzation 
probability Pn{El) of the network over S in terms of the 
orders of growth of the number of L and NL-nodes over S. 

Theorem 1 (dense networks). Let S be a bounded disk of 
radius R belonging to K^. Assume that two sets of nodes with 
statistically independent Poisson point processes with densities 
PL and pjvl are deployed over 5 C 3fi^. Let Nl and Nml 
be the number of L-nodes and NL-nodes, respectively, falling 
in S, and assume that Nl and Nml asymptotically grow as 
the functions fL{n) and fNL{n), where n is an asymptotic 
growth parameter. 

The network of NL-nodes gets a.a.s. localized, i.e., 

lim Pn{El) = 1 

n — >-oo 

for any fhi'n) and fNL{n) such that 

lim /7VLW/l(n)e-^/-(") =0 

n — *oo 

whereby 7 is an appropriate real constant greater than zero. 

Proof. Consider (El) in ( fTSl l along with the relation ( fTTl i. 
and the following inequalities [27]: 

(1 + x)" < e"^, Vxen, x^O (27) 
1-xy < {1- x)y, < X <1 <y (28) 

Based on the previous two relations, Pn (El) in ( fTSl l can be 
bounded as follows: 



1 - X{Xnl) ■ Nnl < Pn (El) < e 



•A' A 



(29) 



where, Nnl > 1 and X{Xnl) < 1 by definition. Equ. 
will be used for demonstrating the three claims of the theorem. 

It suffices to demonstrate that as n 00, X{Xi\fL)-NNL — > 
so that Pn (El) 1, i.e., the network of NL-nodes over S 
gets localized w.h.p. 

Let us rewrite X{Xnl) in an appropriate form for succes- 
sive developments. Upon setting 



(30) 



7 = [^) 
Nl = PlttR^ 

X{Xnl) ■ Nnl can be rewritten as follows: 

XiXNL) ■ Nnl = iVjvLe-^^^ [l + jNl + ^I'^Nl] 
= c- NnlNIc-'^^^ 

(31) 

whereby c = + ^ + ^] . 

From dSTJ, it is straightforward to demonstrate that for any 
fL{n) and fNL{n) such that 

lim NnlNIc-'^^'- = lim fNL{n)fl{n)e-^^'-^'''> =0 

n — >oo n — >oo 

the network of NL-nodes over S gets localized a.a.s. 

□ 

The previous theorem is the starting point for identifying 
appropriate orders of growth of both L and NL-nodes 
guaranteeing asymptotically almost sure localization. In this 



respect, we note the following corollary. 

Corollary (dense networks). Under the scenario described in 
Theorem 1, as n ^ 00 the following holds: 

1) Suppose Nnl - fNL{n) - q ■ n^-^ with ^ G [0, 1) 
and Nl ^ fL{n) ~ p ■ ln(n), with p and q two suitable 
constants strictly greater than zero. 

Then, the network of NL-nodes over S gets localized 
w.h.p. as n ^ c» provided that 

2) Suppose Nl ^ fL{n) ^ In {fNL{n)). 

Then, the network of NL-nodes over S gets localized 
w.h.p. as n ^ 00 provided that 



> 1 



R 



3) Suppose Nl ^ fL{n) ^ n and Nnl ^ fNL{n) ^ n* 
with t > as n ^ 00. Then, the network of NL-nodes 
over S gets localized w.h.p. as n ^ 00. 

4) As a consequence of the previous point, suppose Nnl — 
fNL{n) ~ 0(1), that is, Nnl is a bounded sequence. 
Then, the network of NL-nodes over S gets localized 
w.h.p. as n ^ 00 provided that Nl ~ ^ uj{n) 
with u!{n) ^ 00 no matter how slowly uj{n) grows. 

Proof. As far as claim 1) of the corollary is concerned, it 
suffices to demonstrate that as n 00, X{Xnl) ■ Nnl 
for Nnl ^ qn^^^ with ^ e [0, 1) and Nl ^ pln{n) with p 
and q two suitable constants strictly greater than zero. 

If Nl ~ pln(n) + o(ln(?i)) with p a suitable constant/) > 0, 
it follows that, 

X{Xnl)-Nnl - c-pMn2(n)7VA,i.e-^P'"(") 
= c-p^\n'^{n)NNL-n-''P 



(32) 



In the case Nnl ^ 
have: 

X{Xnl) ■ Nnl 

When n 00, X{Xnl) ■ Nnl 
holds: 

1 - ^ - 7P < 



qn^ ^ with ^ G [0, 1), for n ^ 00 we 



c-q-p^\n{n)-n^-^-^P (33) 
if the following relation 



since we have [27], 



lim 

X — >oo 



MxT 



rP 



0, V > 



By substituting the definition of 7 in the previous relation, 
after some algebra the following threshold follows: 

P>Po= (34) 



Claim 2) follows from observing that for Nl ^ fL{n) 
In {fNL{n)), (ET) can be rewritten as 

\l-7 



{fNL{n)Y~^\n'{jNL{n)) (35) 



X{Xnl) ■ Nnl 
As n — > 00, it is 

lim {fNL{n)f'nn^{fNL{n)) = 



10 




Fig. 1 1 . Behavior of tlie localization probability Pjv {El ) as a function of 
the constant p in A'^^ p ■ log{n) for unboundedly increasing values of n. 
Transmission scenario is compliant with the following parameters; ^ = 0.51 



(Nnl 



9 dB, n„ 



4, I3th = 30 dB, R = 60 m, and 



; 5.62 m. With this setup, the threshold po = 36.72. 



provided that 1 — 7 < 0, from which 



7 ^ 



R 



> 1 



Claim 3) follows from observing that for Nl 
Nnl ~ the following holds; 



X{Xnl) ■ Nml 



0, n 



n and 



(36) 



no matter what the order t of growth of the number of 
NL-nodes. So, asymptotically, the network of NL-nodes gets 
always localized w.h.p. under these conditions. 

Finally, claim 4) follows from the proof of claim 1) upon 
considering = 1 in dSTl i. Note that based on the proof 
of claim 1), ^ = 1 signifies the fact that Nnl = 0{1), 
i.e., Nnl is a bounded sequence, and that X{Xnl)Nnl ^ 



Nle~'^^^ for any Nl ^ w(n) 



cxD as n 



Since inequality ( 1341 ) in Claim 1) is the most important result 
of this corollary, some considerations are in order The basic 
meaning of this result is as follows; in a bounded circular 
region 5 C with area ttR^ with R ^ dmax, the network 
of randomly deployed NL-nodes gets asymptotically localized 
even though the number of L-nodes grows only logarithmically 
(i.e., with an order of growth smaller than that of the NL- 
nodes) provided that the constant p is above the threshold 
Po- This result is fundamental from a point of view of 
network topology, since it assures us that a number of L-nodes 
which grows only logarithmically suffice for assuring network 
localization, provided that p > pq, even though the number 
of NL-nodes grows faster than logarithmically. It is worth 
noting that these results also hold for random geometric graphs 
(RGG); in a transmission scenario typical of RGGs, whereby 
any NL-node can communicate with any other L-node within 
the distance r = dmax, we have cTs = (=J> r] = 0), and the 
threshold becomes: 



PO,RGG 



R 



Borrowing the terminology used in the context of random 
graph theory [33], claim 1) of the previous corollary states that 
the function Nl ^ po ln(n) + o(ln(n)) is a threshold for the 
localization problem at hand. Any function Nl ^ o{po ln(n)) 
allows network localization asymptotically w.h.p. 

As a reference example. Fig. [TT] shows the behavior of 
the localization probability Pn{El) as a function of p for 
unboundedly increasing values of n in the transmitting sce- 
nario summarized in the figure caption. Note that, for p < 
Pq = 36.72, Pn{El) is always zero, while Pm{El) becomes 
instantaneously unitary so long as p = po while n ^ 00. 

Finally, notice that such a threshold does not hold for 
single NL-node localization probability. In other words, upon 
considering the probability P{El) in ( fTSl l for single NL-node 
probability, it is simple to observe that any randomly chosen 
NL-node over a bounded domain S gets localized w.h.p. for 
Nl ~ Lo{n), whatever the behavior of the function Lo{n), 
provided that Loin) -^00 as n ^ 00. 

The results obtained for dense networks state conditions for 
a.a.s. localization of a network of NL-nodes over a bounded 
circular domain for a variety of orders of growth of the number 
of NL-nodes deployed. 

Let us now look at the problem from a different perspective. 
In other words, we look at the problem by considering constant 
L-node density while we let the size of the domain S to grow 
in such a way that pL — — 0{1). Such a result is typical 
of non-dense networks. In this respect, it is useful to evaluate 
the minimum dmax above which the network of NL-nodes 
gets localized a.a.s. 

Theorem 2 (unbounded domains, constant densities). Let S 

be a disk of radius R belonging to 5R^. Assume that two sets 
of nodes with statistically independent Poisson point processes 
with densities pL and p^L are deployed over S* C 5R^. Let Nl 
and Nnl be, respectively, the number of L-nodes and NL- 
nodes falling in S, and consider any asymptotically increasing 
function oj{n), such that uj{n) ^ 00 as n — > 00, and assume 
that Nnl o(w-2(„)e+'^(")). 

Moreover, assume that, as i? ^ 00, the L-node density 
satisfies the following relation: 



PL 



0(1) 



(37) 



Then, as n 00 in such a way that ( [37] i holds, the network 
of NL-nodes gets a.a.s. localized if. 



-LL!{n) 



(38) 



Proof. The proof follows an outline similar to the one of 
the previous theorem. Consider Pn{El) in (fTSb along with 
its bound in (l29b . As before, the objective is to show that 
asymptotically, the transmission range d,nax between each pair 
of L-NL-nodes should grow at least as specified in (l38l) in 
order for Pn{El) 1 as n 00. 
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Given TVtvl, X{Xnl)Nnl can be rewritten as follows: 
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X{Xnl)Nnl = Nnlc-p^^^ 
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With this setup and given (|29l l, it suffices to show that 
X{\nl)Nnl when dmax grows as stated in ( |38] |. 

Upon substituting dmoj: given in dJSl l in (|39] l, the following 
relation follows: 

^(AjvL)iVjVL = iVjVL • e-"(") [1 + oj(n) + 

which goes to zero so long as Lo{n) ^ 00 as n cx) for 
any Nnl — o{u}^^{n)e^^^^'>), guaranteeing that the network 
of NL-nodes gets localized w.h.p. 



The result stated in this theorem is reminiscent of percolation 
theory. In other words, when the deployment region S tends 
to become the entire plane (i.e., R — > 00) in such a way 
that PL is a finite and constant value, the entire network of 
NL-nodes becomes a giant localized component so long as 
the transmission range dmax takes on the values expressed by 
provided that Nnl = o(cj"2(„)gi^(n))^ 

As an example, if uj{n) ^ ln(n), and 
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Iv? n 



the network with an ever-increasing size gets asymptotically 
localized so far as d^^^ grows at least as d^ax ^ 

Notice that, since in practice no real device can support an 
ever-increasing communication range dmax, as the network 
domain increases in size, in the limit there is always a non- 
zero probability that some node cannot get localized. 



VII. Conclusions 

The aim of this paper has been manyfold. Considering a two 
dimensional domain S C di^ over which two sets of nodes 
following statistically independent uniform Poisson point pro- 
cesses with constant densities pL and pjvL are deployed, we 
first derived the probability that a randomly chosen NL-node 
over S gets localized as a function of a variety of system 
level parameters. Then, we investigated the probability that 
the whole network of NL-nodes over S gets localized. The 
transmission scenario assumed is that of shadow fading. 

Furthermore, we presented a theoretical framework for 
deriving both finite case and asymptotic thresholds for the 
probability of localization in connection with both a single 
non-localized node randomly chosen over the investigated do- 
main, and the whole network of non-localized nodes. Finally, 
we investigated the presence of thresholds on the problem 
at hand for unboundedly increasing values of the number of 
deployed nodes over the domain S. 
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Appendix I 

Upon substituting (|9]l in ( fTOl i. and considering r = c? j , 



1 - erf - In 



7? V '^ma 



R-npL / erf -In 



rdr 

rdr 
(40) 



By employing the substitution y — ^ In ^ '' ^ ^ r — 

dmaxC^^ , from which dr = dmax-^e^^dy, the integral ( |40] i 
takes on the following form: 



erf I - In 



rdr^dl^^^^ I eif {y)e^'^ydy 



whereby, -^s = ^ In ( 



Upon using the following [34] 



e'"^erf(6a;)dx = - 



a t (I 

e""^ erf (5a;) - eiP'erf [hx ~ — 



after some algebra, ( |40| ) can be rewritten as follows: 



Eid^'^m = ^R'~np,^crf(^\n 



^d^ 



1 + erf I - In 



R 

} 

'max 

R 



(41) 



Next consider evaluating E{d^^} over 3?^. In the limit R 
oo, i4T[ simplifies to: 



Eid^"^} = lim Eid^'^lR} = pLird 



(42) 



since. 



lim erf(a;) = 1 



